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The Argand diagram is used to display some characteristics of the Riemann Zeta function. The zeros of the Zeta function 
on the complex plane give rise to an infinite sequence of closed loops, all passing through the origin of the diagram. This leads 
to the analogy with the scattering amplitude, and an approximate rule for the location of the zeros. The smooth phase of the 
Zeta function along the line of the zeros is related to the quantum density of states of an inverted oscillator. 

PACS numbers: 05.45.+b, 03.65.Nk, 02.90.+p 



On The Riemann Zeta function C,{s) of the complex vari- 
s = IT -\- it has an infinite number of zeros on the 
"half-line a = \ [I]. These zeros are of great interest to 
^mathematicians from the number theoretic point of view 
^^[2,3], and to physicists interested in quantum chaos and 
the periodic orbit theory [4-6]. Along this line, as a func- 



jtion of t, every time C(s) changes sign, a discontinuous 
jump by TT in the phase angle is introduced. Otherwise 
the phase angle is a smooth function of t. The smooth 
^~^art of the phase angle itself is very interesting, since it 
f~younts the number of zeros on the half-line fairly accu- 
V^rately. In this letter, we bring out some of the interesting 
^Dproperties of (^(s) by displaying it on an Argand diagram, 
where 3'iC(s) (along the x-axis), is plotted against 9C(s) 
^i,^(y-axis). For a = 1/2, the plot yields an infinite sequence 
Gof closed loops, one for every zero of the Zeta function, all 
^! &oing through the origin. This character changes dras- 
I tically as one moves away from the a = 1/2-line. The 
Osimilarity in the Argand diagrams for C(l/2 + H) and the 
J^resonant quantum scattering amplitude is pointed out, 
""^and this analogy, although flawed, leads directly to an ap- 
proximate quantization condition for the location of the 
zeros [7,8]. We also demonstrate that for a = 1/2, the 
smooth part of the phase angle of C(s) is closely related 
to the phase shift of an inverted half harmonic oscillator. 

We start by giving some standard results for the Rie- 
mann Zeta function C(s). Using the fundamental func- 
tional relationship between C(s) and C(l — s), it is easy 
to show that 



where Z{t) is real, and 9{t) is the phase angle, with 
the convention that (^(O) = 0. Using Eqs.(l) and (2), it 
follows that 



exp(2i6'(t)) = exp(— itluTr) 
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The phase 9, as defined above, is smooth in the sense 
that it does not include the jumps by tt due to the zeros 
oi Z{t). Nevertheless, the number of zeros between and 
t on the (7=1/2 line is counted fairly accurately by 0{t), 
as will become clear from the Argand diagram. Note that 



T^=-2^^^'^+2^^[^^^(4+2)-^^^(4-2^ 



The density of zeros is given by 



TT dt 27r 
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where the digamma function is defined as ^(2) = 
r (z)/r(z). From the above, the asymptotic expression 
for 6(t) may be obtained immediately by making asymp- 
totic expansion of the F functions. We denote this by 
6(t) , and it is given by 



TT ^ ' ^2tt' ^2tt' ^2tt' 8 487rt 
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where F(z) denotes the gamma function of the argument 
z. We may further write 



a- + it) = Z{t)exp{-ie{t)), 



(2) 
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behaviour of the phase Bi^) away from the (7 = 1/2 line. 
This is also apparent in Figs. 2(a) and 2(&), where the 
phase angle Bi^) is plotted as a function of t on and off 
the 1/2-axis. The abrupt change in the behaviour of the 
phase angle as shown in Fig. 2, is a consequence of 

the fact that for a ^ 1/2, C(s) and C(l — s) are not com- 
plex conjugate of each other. Therefore Eq.(3) no longer 
holds in such a situation. The "chaotic" characteristic 
of the phase angle off the 1/2— axis, is, of course, well- 
known. 
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FIG. 1. The Argand diagrams for ^(cr -|- it) for increasing 
t and fixed a. (a) C(l/2 + it) in the range t = 9 — 50. The 
lower limit for t is chosen so as not to miss the Gram point 
flanking the lowest zero at t = 14.13. (b) C(l/2 + it) for 
t = 280 — 300, showing the two loops without the correspond- 
ing Gram points, (c) C(0-6 -|- it) for t = 9 — 50 to show the 
defocussing at the origin, (d) ^(1 -|- it) for t = 9 — 50. Note 
the pronounced shift of the diagram away from the origin in 
this case. 
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To bring out some characteristics of the function 
C(l/2-|-it), we plot its Argand diagram in Fig. 1(a) in the 
range t = 9 to t = 50. This shows a sequence of closed 
loops, one for every zero of the Zeta function. At a zero 
of C(l/2 + H)-, both its real and imaginary parts are zero 
at the same value oft, and therefore every loop converges 
at the origin. The intercepts on the real axis are the so- 
called "Gram points" , where only the imaginary part of 
C(s) is zero, due to the phase angle 9{t) = nvr. With 
infrequent exceptions, there is one Gram point between 
two consecutive zeros of the Zeta function. The first two 
exceptions to this rule occur for the 126th and the 134th 
zeros at t = 282.455 and 295.584 respectively [2]. The 
Argand diagram in Fig.l(&) clearly shows that for these 
cases, the loop structure still persists, even though the 
Gram points are missing. In Figs. 1(c) and l{d), Argand 
diagrams are drawn away from the 1/2-axis, for a = 0.6 
and (7 = 1 respectively. These clearly show the defo- 
cussing at the origin due to the absence of the zeros in 
the Zeta function. Moreover, the number of intercepts 
along the real-axis in the Argand diagrams now show a 
large increase compared to the (7=1/2 case, whereas 
the intercepts on the imaginary axis are few or nonexis- 
tant. This is a reflection of the dramatic change in the 



FIG. 2. The exact phase angle of C{a + it) is plotted as a 
function of t = — 50 for (a) (7 = 1/2 and (b) a = 0.6. See 
text for discussion. 

Indeed, the phase angle for s = 1 -|- it is well studied, 
and is related to the quantum scattering phase shift of 
a particle on a surface of constant negative curvature 
[13-15]. 

Finally, in Figs. 3(a) and 3(&), the Argand diagrams 
of the Zeta function are drawn for a much larger range 
of t, from 1 to 500, on and off the 1/2— axis. Note that 
the scale for (7 = 1 is expanded compared to that for a = 
1/2. Borrowing from the terminology of the motion of a 
particle in phase space, it is as if there is an "attractor" 
at the origin for (7=1/2 (Fig. 3(a)), which is absent 
from the more disorderly tracks of Fig. 3(&), which is 
drawn along the (7=1 line. The latter figure also shows 
that the real part of the Zeta function is always positive 
for (7=1 for this entire range oft. 
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f,{k) = {tj, exp{2i6,) -l)/2ik. 



(7) 
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FIG. 3. The Argand diagrams of the Riemann Zeta func- 
tion for the wider range of t— values, from t = 9 — 500. (a) 
C(l/2 + it), and (b) C(l + 

The loop structure of the Zeta function at a = 1/2, 
with some near-circular shapes, is reminiscent of the Ar- 
gand plots for the scattering amplitudes of different par- 
tial waves in the analysis of resonances, for example in 
pion-nucleon scattering [9]. Consider the partial wave 
amplitude fi{k), defined in terms of the partial wave 
phase-shift 6i[k) and the inelaticity parameter rji^k), 



Here / refers to the angular momentum, and k the wave 
number. Note that 9/; (A;) is never negative, since the 
inelasticity parameter rji is always less than one. One 
generally plots an Argand diagram with 2kQ fi[k) along 
the J/— axis and 2¥iR fi{k) along the x-axis for various val- 
ues of k. For the case of no inelasticity (ry; = 1) and a 
single resonance, the Argand diagram is a perfect circle 
with unit radius, with the center on the imaginary axis 
at 1. By comparing this with Fig. 1(a) at cr = 1/2 , we see 
that the real and the imaginary parts are interchanged 
in the latter, but otherwise there is a strong similarity, 
with many of the loops having inelasticity. This analogy 
is flawed, however, since ■KC(l/2 + it) does become neg- 
ative in small islands of t. Nevertheless, if these islands 
are ignored, then the phase shift bi may be identified with 
the phase angle 6'-|-7r/2, with each closed loop in Fig. 1(a) 
being regarded as in isolated resonance . In this approx- 
imation, the Gram points occur as before for sin 6* = 0, 
while the zeros of C(l/2 + H) are given by the condition 

cos6' = 0, 6* = (m + l/2)7r, m=l,2,.... (8) 

This condition for the location of the zeros was also ob- 
tained by Berry [7] from the first term in his approxi- 
mate formula. Eq.(8) has roots that yield the zeros on 
the 1/2— axis with an error of at most 3 percent. 

As mentioned earlier, for cr = 1, the phase 9{t) is inti- 
mately connected to the quantum scattering of a particle 
on a saddle-like surface. On the (7=1/2 line, our anal- 
ogy with the scattering amplitude also suggests that the 
phase angle 9{t) is related to a scattering phase shift. 
We now demonstrate that the scattering of a nonrela- 
tivistic particle by an inverted harmonic oscillator with 
a hard wall at the origin generates a phase shift that is 
closely related to 0{t). Indeed, we show that the quan- 
tum density of states for this problem is essentially the 
same as Eq.(5) for the density of the zeros. Consider the 
Schrodinger equation for a; > 0, 

-- — - -mw^a;^* = (9) 
2mdx^ 2 ' ^ ' 

and impose the boundary condition that the wave func- 
tion $ vanishes at the origin. Putting = y, ^ = y~i(j>, 
it becomes 

£ 1(1 + 1) kt ^ 

dy^ y^ y 

In the above equation, 

1 , E , ^ ^ , 

This is effectively a three dimensional Schrodinger equa- 
tion for a repulsive Coulomb potential in the variable y. 
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To obtain the phase shift, we write the asymptotic solu- 
tion of the above equation [16] as 

^{y)^sm{ky-Un{2ky)-'-^ + 7j,), (12) 

where rji is the phase shift with respect to the distorted 
Coulomb wave, given by arg r(/+ 1 + it/2). For our one- 
dimensional problem, only / = — 1/4 is relevant. For this 
case, omitting the subscript /, the phase shift rj is 



„/3 it^ 
V = argr(- + -). 



Using the identity [10] 
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T(- + iy)r( iy) = — . 
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the number of quantum states n(t), between and t, is 
then given by 

nW=^=^ + ^9 [inr(- + .-)-inr(- -.-)]. 

(15) 

In the above equation, C is a smooth function given by 



C 



TT 

2 



tan ^(cosechTrt ). 



(16) 



Note from above that the expression for ry(t) is not quite 
identical to 9{t) as defined in Eq.(4). However, their 
derivatives, the quantum density of states, only differ by 
a constant and an exponentially small term. It should 
also be pointed out that even if we had started with a 
full inverted harmonic oscillator (rather than the half- 
oscillator), the same conclusion would be reached, even 
though there may arise some nonuniqueness in the choice 
of the boundary condition. The inverted harmonic oscil- 
lator problem was studied by a number of authors in 
the past [11,12] in relation to time-delay, and by others 
[17-20] in connection to string theory. No connection, 
however, was made to the phase of C(l/2 + H)- 

In summary, we recapitulate the main points made in 
this paper. Traditionally, the behaviour of the zeros of 
the Riemann Zeta function on the 1/2— axis is associated 
with the bound state problem of a quantum Hamiltonian. 
We, on the other hand, focus on the scattering problem. 
We use the Argand diagram construction for the Rie- 
mann Zeta function to visualize how the smooth phase 
9{t) acts as a counter for the zeros of Z{t). This leads 
to the analogy with the scattering amplitude, and the 
approximate condition ( Eq.(8) ) for the location of the 
zeros. We also find a connection between the phase angle 
of the Zeta function on the (7=1/2 line and the quantum 
phase shift of a potential. It is perhaps not surprising 
that this potential is as simple as the inverted harmonic 
oscillator, and not something more complicated, as is the 



case for the chaotic phase for (7=1 [13-15]. Finally, we 
suggest that the Argand diagram construction may also 
be useful for the Selberg Zeta function, and shed light on 
the quantisation condition for quantum chaos [8]. 
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